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Abstract 

We review the applications of the integral over anticommuting Grassmann variables (non- 
quantum fermionic fields) to the analytic solutions and the field-theoretical formulations 
for the 2D Ising models. The 2D Ising model partition function Q is presentable as the 
fermionic Gaussian integral. The use of the spin-polynomial interpretation of the 2D Ising 
problem is stressed, in particular. Starting with the spin-polynomial interpretation of the 
local Boltzmann weights, the Gaussian integral for Q appears in the universal form for a 
variety of lattices, including the standard rectangular, triangular, and hexagonal lattices, 
and with the minimal number of fermionic variables (two per site). The analytic solutions 
for the correspondent 2D Ising models then follow by passing to the momentum space on 
a lattice. The symmetries and the question on the location of critical point have an inter- 
esting interpretation within this spin-polynomial formulation of the problem. From the 
exact lattice theory we then pass to the continuum- limit field-theoretical interpretation of 
the 2D Ising models. The continuum theory captures all relevant features of the original 
models near T^. The continuum limit corresponds to the low-momentum sector of the ex- 
act theory responsible for the critical-point singularities and the large- distance behaviour 
of correlations. The resulting field theory is the massive two-component Majorana theory, 
with mass vanishing at Tc- By doubling of fermions in the Majorana representation, we 
obtain as well the 2D Dirac field theory of charged fermions for 2D Ising models. The dif- 
ferences between particular 2D Ising lattices are merely adsorbed, in the field-theoretical 
formulation, in the definition of the effective mass. 
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Abstract 

We review some aspects of the anticommuting integrals as applied to the analytic 
solutions and the field-theoretical formulations for the 2D Ising models. We stress, 
in particular, the use of the spin-polynomial interpretation of the 2D Ising problem. 

1. Introduction. A remarkable feature of the two-dimensional Ising model, first es- 
tablished within the transfer-matrix and combinatorial considerations, is that it can be 
reformulated as a free fermion field theory. It was recognized later on that the natural 
tool to analyze the 2D Ising models is the notion of the integral over Grassmann variables 
(nonquantum fermionic fields) introduced by Berezin [1-8]. A simple interpretation of the 
2D Ising model based on the integration over anticommuting Grassmann variables and 
the mirror-ordered factorization principle for the density matrix was introduced in [7,8]. 
The partition function is expressible here as a fermionic Gaussian integral even for the 
most general inhomogcncous distribution of the coupling parameters [7]. This also gives 
a few line derivation of the Onsager result [7]. The Gaussian fermionic representations 
has been recently constructed as well for the inhomogeneous dimer problems [9]. The 
method [7,8] do not involve traditional transfer-matrix or combinatorial considerations. 
Schematically, we have: 

Q = Sp Q(a) ^ Sp Q (a\a) ^ Sp Q (a) = Q . /-in 

(a) (a I a) (a) ^"^^ 

We start here with the original spin partition function Q, and then we introduce, in 
a special way, a set of new anticommuting Grassmann variables (a), thus passing to a 
mixed {a \ a) representation. Eliminating spin variables in this mixed representation, we 
obtain purely fermionic expression (fermionic Gaussian integral) for the same partition 
function Q. Grassmann variables are introduced, at the first stages, by factorization of 
the local weights. The factorization ideas resemble, in a sense, the idea of insertion of 
Dirac's unity, E ja )( a| = 1, in transformations in quantum mechanics [10]. An important 
ingredient of the method is also the mirror-ordering procedure for the arising noncommut- 
ing Grassmann factors [7,8]. The original variant of the method implies two Grassmann 
variables per bond in the resulting fermionic integral, since we start here with the fac- 
torization of the bond Boltzmann weights [7] . An important modification was introduced 
in [8], where we start with the factorization of the cell weights presented by three-spin 
polynomials. This enables us to obtain the fermionic representation for Q with only two 
fermionic variables per site which provides essential simplifications in the analysis [8]. 
The symmetries and the criticality conditions have interesting spin-polynomial interpre- 
tation [8]. In this report we comment on the analytic results obtained so far within this 
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spin-polynomial fermionic approach and derive the effective field-theoretical (continuum- 
limit) formulations for the 2D Ising models valid near T^. The resulting field theory is 
the massive free-fermion Majorana theory. The form of the action is universal, particular 
Ising lattices, like rectangular, triangular and hexagonal lattices, are specified merely by 
the definition of mass. By doubling of fermions in the Majorana representation, we obtain 
as well the 2D Dirac field interpretation for 2D Ising models near Tc. 

2. Free-fermion representation and analytic results (lattice case). 

The spin-polynomial interpretation of the 2D Ising models can be illustrated by an 
example of the triangular lattice [8]. Let us consider the triangular Ising lattice as a 
rectangular one with the additional diagonal interaction introduced in each lattice cell. 
Let mn be the lattice sites on the corresponding rectangular net, m,n = 1,2,3,...,L, 
where L is lattice length, N — L'^ is the number of lattice sites. At final stages we assume 

— > oo. The hamiltonian of the triangular 2D Ising model is [8]: 

-P H (a) = J2 [ bl<ymn<^m+ln + b2Crm+ln<^m+ln+l + h<^mn<^m+ln+l ] , P = l/kT , (2) 
mn 

where amn — il are the Ising spins, ba — Ja/kT are dimensionless coupling constants, Ja 
are the exchange energies, kT is temperature in energy units. In what follows we assume 
the purely ferromagnetic case, ba > 0. Respectively, the partition function and the free 
energy per site are given as follows: 

^ = E(.)e-/^^(-), -/3/z= Jhn ^InZ, (3) 

where the sum is taken over 2-^ spin configurations provided by amn = ±1 at each site. 
Noting the identity for the typical bond weight e^'^'^ = cosh b + sinh b-aa', which readily 
follows from (aa')^ = we obtain: 

Z = (2 cosh bi cosh 62 cosh 63) ^ <3 , = ^ 00 , (4) 

where Q is the reduced partition function [8]: 

Q — Sp { n ( 1 + tlCmn(^m+ln) (1 + hCm+ln^m+ln+l) (1 + t3<^mnCm+ln+l) } , ( r.\ 

(o-) mn 

where t^ — tanh ba, and Sp(^a) — 2~^S(o.) is the normalized spin averaging. 

The spin-polynomial interpretation arises if we multiply the three bond weights form- 
ing a triangular cell in Q. Noting that aia2 ■ j, = (Jia^, etc, we obtain [8]: 

Q — Sp Yl {oiQ + aiamnO'm+ln + Ol2 0'm+lnO'm+ln+l + Ol3amnO'm+ln+l) , (Q) 
(u) mn V / 

where for the triangular lattice we have ckq = l-\-tit2t3 , cci = ^1+^2^3 , 0:2 = ^2+^1^3 , 0:3 = 
^3 + tit2 , the case of the rectangular lattice follows with ^3 = 0. As is discussed in [8] 
the hexagonal and some other lattices are covered by (6) as well, with the correspondent 
specification of the a parameters. In what follows we consider the three-spin polynomial 
partition function Q, assuming that a^, ai, a2, as are some numerical parameters. 
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The partition function (6) can be transformed into the fermionic Gaussian integral as 
follows [81: 



ran mn ( / ) 



where c^„, c*^^ is a set of purely anticommuting Grassmann variables. The free-fermion 
representation for Q given in (7) is exact and completely equivalent to the original spin 
representation (6) up to the boundary effects negligible in the limit of infinite lattice 
— s> oo. [We do not comment here neither on the properties of the Grassmann variables, 
nor on the method of the derivation of this representation. The derivation is simple and 
takes one page of transformations, neglecting preliminary preparations, the details are to 
be seen in [8] . Also see the comments on the fermionic Gaussian integrals in [9] ] . 

The explicit evaluation of the integral (7) can be performed by passing to the momen- 
tum space. For the squared partition function we find [8]: 



L-l L-1 



<9^=n n [ ("o + + «i + «3) - 2 (Qio«i - «2«3) cos ^ - 

p=o <}=o (8j 

— 2 {aQa2 — a^ag) cos 2|2 _ 2 {a^a^ — a^a2) cos ^^^^'^^ ] . 
Respectively, for the free energy we obtain: 

lim^^^(^lng) = 

27r 27r 

= i / / f 1^ 1" [ («o + + «i + «i) - 2 (aoai - a^a^) cosp - (9) 

— 2 {aQa2 — Oiia^) cosg — 2 {a^a^ — q;jQ;2) cos(p + q)] ■ 

It is easy to check that the free-energy form (9) yields the known exact expressions for 
the free energies of the rectangular, triangular, and hexagonal Ising lattices, with the 
correspondent specifications of the parameters . 

3. The symmetries and the critical point. 

The symmetries provided by the above exact solution and the question on the location 
of the critical point have an interesting interpretation in the spin-polynomial language [8] . 
Introducing the notation {ai,a2,a3)mn ((Jmn, c^^+in, o-m+in+i), the cell weight in the 
density matrix in (6) is given as the three-spin polynomial: 

-Pi23 (c) — ao + ai (7i(72 -I- a2 (72(73 -|- 0:3 (71(73 . (10) 

It is also useful to introduce the associated three-spin polynomial: 

-^123 (cr) ^ (Xq- ai (7i(72 - a2 (72(73 - (7i(73 , (11) 

and we note the following interesting identity [8] : 

(Fi23((7)) ^ = (C^o - "1 0"l<72 - <^2 (^20"3 ^ "3 0"lO"3 ) ^ = 

= {al + al + al + aj) - 2 (ao"! - '3;2'3;3) cric^2 - (12) 
— 2 {aoa2 — aias) (72(73 — 2 (q;oQ;3 — aia2) (71(73 . 
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It is seen that the combinations of the a-parameters occurring in (^123)^ are exactly the 
same as in the momentum modes Q'^{p\q) in (9): 

Q"^ {p\q) = {al + al + al + af) — 2 {a^a^ - Q;2«3) cosp- . . 

— 2 (q;o'^2 ~ CKj^ctg) cos q — 2 (ctQCCg — aj^ag) cos {p + q) , 

where < p, g < 27r, the limit ^ 00 is already assumed. 

It appears that the following combinations of the a-parameters play important role in 
discussing the symmetries and the critical point [8] : 

Q!o = I (ao + «! + ^2 + eta) , 
Q!j[ = I (cto + ai — a2 — eta) , 
^2 — \ (<^o — ai + 0:2 — 0:3) , 
^3 = I (tto - ai - 0:2 + 0:3) , 

The parameters a* and a are in fact the eigenvalues of the polynomials \P122, and |Fi23, 
respectively. By the "eigenvalues" we mean the four numbers which takes the polynomial 
as the spin variables run over their permissible values ±1. We note also the following 
identity: 

0.^0.-10.20.^ — 0^01020"^ — q;oQ;iQ;2Q;3 . (15) 

There are some evident symmetries in the solution. For instance, the free energy (9) is 
a symmetric function with respect to arbitrary permutations of ao, ai, 02, parameters. 
We can as well change the signs of any two of them, with —^fg unchanged. There is also 
a less evident hidden symmetry in the solution, corresponding to the Kramers- Wannier 
duality in the case of the standard rectangular lattice. Namely, the partition function 
Q{a} is invariant under the transformation Oj o*. This symmetry in fact holds 
already for the parameters of the separable fermionic modes (13), and can be proved 
making use of (12), see [8]. 

In order to establish the possible critical points, we have to look for zeroes of {p\q) 
momentum modes. As it can be guessed already from the analogy between (12) and (13), 
there are four exceptional modes with {p\q) = (0 | 0), (0 | tt), (tt | 0), (tt | tt). For these 
modes we have, respectively, {p\q) — (2ao)^ , (2q;i)^ , (20:2)^ , (20:3)^ . We remember 
that the parameters and hence Oj are some functions of temperature. Thus, if at some 
temperature one of the above momentum modes vanishes, we fall at the point of phase 
transition. It can be shown that all other modes Q^{p\q) are definitely positive at all 
temperatures, there are no other critical points for physical values of the ct-parameters. 
The possible criticality conditions can be combined into one equation: 

Oioaia20i2 = 0. (16) 

It can be shown also that if all bond interactions are ferromagnetic then ai, 0:2, 0:3 never 
become zero, and the critical point can only be associated with ao = 0, or cto — cti — 0:2 — 
03 = 0. The criticality conditions with = 0, A; = 1, 2, 3, can only be realized when the 
antiferromagnetic interactions are involved in the hamiltonian. 



"0 = 

Oil = 

02 = 

Or, — 



ao — «! — tt2 ~ <^3j 



2 (cto — «! + q;2 + 0^3) , 
I (CKO + CCi — 0:2 + o-^ , 
\ (cKo + CCi + 0:2 — 0:3) . 



(14) 
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If the critical point is associated with aj — 0, then the singular part of the free energy 
(9) near is given as follows [8]: 



(2aj)^ ^ const 



-(^Iq I singular " — / , ^ = 77^^~yi ' i^'^) 



which implies that near Tg in order of magnitude —Pfq ~ Inr^, where r ~ \T — Tc\. 
The specific heat have thus the log-type singularity, C ~ | lnr| as T Tc (Onsager, 
1944). 

It is seen that the eigenvalues (14) play important role, but it is not so clear, in fact, 
how the role of the polynomial F123 can be understood in a less formal way, at the level of 
the original spin-variable formulation of the problem, prior to the analytic solution. Even 
more amusing grounds to be interesting in this game with the parameters a^, a^, a* 
follow from the expression for the spontaneous magnetization below T^. The 8-th power 
of the spontaneous magnetization M —< Umn > for model (6) is given by the following 
very simple expression [8]: 

^8 ^ r-^-. «o«i«2a3 ^ _ aoala*2a*3 ^^^^ 

This expression for holds true when the right hand side varies between and 1, and 
— otherwise. The well known expressions for the spontaneous magnetizations of the 
rectangular, triangular, and hexagonal lattices follow easily from (18) as particular cases. 
Prom (18) we find M ~ rs as r ~ | T — Tc | — 0, with the universal value of the critical 
index /9 = 1/8 for the magnetization at the critical isobar. What are the hidden reasons 
for such simple expression for M^, is yet unknown. 

4- Majorana field theory for the 2D Ising models. 

The fermionic integral for Q is a suitable starting point to formulate the continuum 
field theory near Tc. We write once again the exact lattice action from (7) as follows: 

S{c) =J:[ (ao - ai - Q;2 - as) C^n^mn + «1 ^mn {Cmn " C-ln) + 

mn 

+ Oi2 Cmn (c^n ~ ^mn-l) ^rnn (c^n ~ ^m-ln-l) ^ ' 

with c^^ = c^^ = . Let us define lattice derivatives dm Xmn = Xmn — Xm-in , dn Xmn = 

Xmn - Xmn-1 , with Xmn " Xm-ln-1 = dm Xmn + 9„ Xmn " dmOn Xmn ■ Introducing the UeW 

notation for the fermionic fields, , ipmn , i^mn , from (19) we obtain: 

mn _ _ _ 

mn ~l" '4'mn ^n 'ipmn ] ) (20) 

Xi — ai + eta, X2 — Oi2 + (X3 , m — ao — ai — a2 — as — 2 ao . 

It is easy to recognize in this still exact lattice action the typical field-theoretical like 
structure with mass term and kinetic part. Remember that we assume the ferromagnetic 
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case with the criticahty condition 2ao = . The parameter m = cto — cti — q;2 — cts = 2 
play the role of mass in the field-theoretical interpretation, m ~ \T — Tc\ — > as T Tc , 
the ordered phase corresponds to ttt, < . 

Taking in (20) the formal limit to the continuum euclidean space, with (mn) — > 
(a;i|a;2), i/Jmn i>{x) = ijj {xi\x2} and dm ^ di ^ d/dxi, 9„ ^ ^2 = d/dx2, and 
neglecting the second-order 8182 kinetic term, we obtain the Majorana like fermionic 
action of the correspondent continuum theory: 

S^j£x [rn'ijj(x)'^(x) + Ai ijjix) di ^(x) + X2 i^{x) 82 i^{x) + 

(21) 

ai '^{x) di ip{x) + a2 V'(x) 82 i}{x) ] , 

where Ai, A2, m arc defined in (20). A nonstandard feature in this action is the presence 
of the "diagonal" kinetic terms like ipdiip, ipd2ip. This terms can be eliminated by a 
suitable linear transformation of the fields ip.ip. Omitting the details, after a suitable 
transformation of the fields and the d'^x space in (21), we come to the euclidean 2D 
Majorana action in canonical form: 

S = J d'^x [irnipi{x)ip2{x) + ipi{x) 8oipi{x) + iIj2{x) 8Qip2{x) ] , 

(22) 

9o = I {di - 182) , 8q^\{8i-\- 182) , 2D Majorana , 



with the rescaled mass: 

a{) — ai — a2 — c 

m — 

1/2 



m = -rj^ . (23) 



(2 ^(00010203)0) 

The 2D Ising model is presented in (22) as a field theory of free massive two-component 
Majorana fermions over the euclidean d'^x space. The Majorana fields V'l ■ ^'^'2 in this 
representation are the linearly transformed fields ip , ip from (21). The axis of the d'^x 
space are also rescaled and rotated as we pass from (21) to (22). Respectively, the mass is 
rescaled according (23), here ( )c means the criticahty condition (cto — ai — a2 — 03)0 = 0. 
The Majorana field theory for the 2D Ising model, for the simplest case of the isotropic 
rectangular lattice, was constructed, by another method, in rcfs. [3,4]. The action (22) is 
a generalization of eq. (129) in [4], also see eqs. (2.11) in the first reference of [3]. The 
generalized action (22) captures the relevant features of the exact lattice 2D1M theories 
for the lattices covered by (6), see also (7), in the low- momentum sector, which is respon- 
sible for the large-distance behavior of the correlation functions and the thermodynamic 
singularities near Tc. In particular, the result (17) can still be recovered from (22). 

5. The Dirac field interpretation of the 2D Ising models. 

We can pass as well to the Dirac field theory of charged fermions by doubling the 

number of fermions in the Majorana representation (22). On formal level this corresponds 
to passing from Pfaffian like Gaussian integral to the determinantal Gaussian integral 
according the identity iPfaff^l^ = det A, see [9]. To this end, we take two identical 
copies S' and 5"' of (22) and consider the combined action dirac — S' + S". In this action 
we introduce the new fermionic fields ipi,ip2, i^i, '02* by substitution: 

V'l = 73(^1 + ^2 ) ,i^i= 7l(V'i - i^i) , 02 = ^(V'2 + 01*) ,i^2= 73(02 - 01*) . (24) 
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In terms of new fields tlie action S = S' + S" becomes: 

-S" dirac ^ Jd^X [im ('ljji(x)ljjl(x) - ljj2(x) ljj^(x)) + 

+ ^i{x) {di - 182) ^*2{x) + ^p2{x) {di + 182) ^l{x) ] , 
where tfie rescaled mass m is given in (23). In matrix notation we find: 



(25) 



^2 



r2 



(26) 



Introducing the 2D cuchdean 7-matrices 71 = Ci, 72 = a2 and 75 = 7172 = icr^, where 
ci,cr2,C3 are standard Pauli matrices, we obtain the same action in the 7-matrix inter- 
pretation: 

S = Jd'^xl'^ {x){ 7n75 + 71 (?i + 72 (?2 ) * ^ {x) ] (euclidean 2D Dirac) , 

(27) 

lulu + lulu = 2 5^^, 7^2 = 72 = 1 , 75 = 7^72, 72 = -1 , 

where ^ = [ipi \ 1^2), ^E'^ = {ip* \ 1P2) arc the two-component charged spinors. We find here 
the euclidean action of the massive free 2D Dirac theory. The unconventional 75 matrix 
in the mass term can be eliminated, if one wants, by a redefinition of the fields, changing 
the sign of a one of the four spinor components. 

We can pass as well to the 2D Minkovsky space. The 2D Minkovsky space is the 
complex plane. Hence we put { a;i | a;2 } { x \ i t} , {xi\ X2) — > ^ (x | it) , di ^ di — 
d/dx, 82^ —ido — —id/dt. From (25) or (26) we obtain: 

-^dirac ^ifd^x^ {x) (as) [m + (^2) + do{ta,) ]^Hx) . (28) 

Introducing the Minkovsky-space 7-matrices 7° = Ui, — —ia2, 7^ = 1^1'^ — ^3; and 
passing to new spinors ^ = ^^7^ , 'l'^ = \1'^, we find: 

S ^ i J d'^x ^ (x) [m + i d] ^^x) (minkovsky 2D Dirac) , 

(29) 

d^ryi^g^^ ryOQ^^^iQ^^ ^ n ^ ^ u ^ ^ 2 Q = dlag ( + | -) . 

We have here the 2D Dirac action in the Minkovsky space interpretation. This is to be 
compared with the action given in eq. (3.17) in [6] for a particular case of the rectangular 
lattice. 

In conclusion, we have discussed the fermionic structure and related aspects of the 
2D Ising models described by the three-spin polynomial partition function (6) and the 
equivalent Gaussian fermionic integral (7). The symmetries and the criticality conditions 
are most naturally expressed in terms of the parameters a, a, a* arising in this spin- 
polynomial interpretation. Also, the spontaneous magnetization is expressible very simply 
in terms of these parameters, see (18), which property is not yet well understood. We 
then derived, directly from the exact lattice action, the continuum-limit Majorana and 
Dirac field theories for the models under discussion, including the standard rectangular. 
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triangular, and hexagonal 2DIM lattices. We find that the lattice parameters are merely 
absorbed, in the continuum limit, in the definition of mass. This confirms the universality 
ideas. Finally, we argue that even the standard rectangular lattice can be most easily 
analyzed, both in the exact lattice theory and in the field-theoretical approximation, 
starting with the spin-polynomial interpretation of the problem. 
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